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The Dirac equation offers a precise analytical description of rela- 
Q ■ tivistic two-particle bound states, when one of the constituent is very 

heavy and radiative corrections are neglected. Looking at the high-Z 
hydrogen-like atom in the infinite momentum frame and treating the 
electron as a "parton", various properties usually attributed to the 
quark distributions in the nucleon are tested, in particular: Bj0rken 
^ , scaling; charge, helicity, transversity and momentum sum rules; ex- 

^ ' istence of the parton sea; Soffer inequality; correlation between spin 

Ph . and transverse momentum (Sivers and Boer- Mulders effects); trans- 

verse displacement of the center-of-charge and its connection with the 
magnetic moment. Deep inelastic experiments with photon or positron 
■ beams at MeV energies, analogous to DIS or Drell-Yan reactions, are 

O ■ considered. 

p 

«N 1 Theoretical frame 

I The Dirac equation enables us to study the relativistic aspects of an 

hydrogen-like atom A of large Z {Za ~ 1, where a = e^/(47r) ~ 
- 1/137). It takes all orders in Za into account but neglects (i) the 

r> i nucleus recoil, (ii) the nuclear spin and (iii) radiatives corrections like 

' the Lamb shift. So it is accurate at least to zero*^ order in a and 

nie/mA- Applying a Lorentz boost, we have an explicit model of 
"doubly relativistic" two-body bound state (relativistic for the inter- 
nal and external motions). In particular, boosting the atom to the 
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"infinite momentum frame" (or looking it on the null-plane t + z = 0) , 
one has a model for the structure functions which appear in deep in- 
elastic scattering on hadrons. In fact, since it neglects nucleus recoil, 
this model is best suited to mesons with one heavy quark. However 
many properties can be generalized to hadrons made of light quarks. 

In analogy with the quark distributions, we introduce the unpo- 
larized electron distributions q{k~^), q{]<.T,k~^) and q{h,k~^) where fe+ 
takes the place of the Bjorken scaling variable, kr is the transverse 
momentum of the electron and the impact parameter b = {x,y) is 
the variable conjugate to k7^. We will also define the corresponding 
polarized distributions like g(b. A;"*", S'^; S"^) where and are the 
polarization vectors of the atom and the electron. We will particularly 
study 

• the differences between q{k'^), the hehcity distribution Aq'(A;+) 
and the transversity distribution Sq{k~^); 

• the sum rules for the vector, axial and tensor charges and for the 

longitudinal momentum; 

• the correlations between S"^, S*^ and b or ky, like the Sivers 
effect; 

• the existence of a non-zero (b) for transverse S"^ and its connec- 
tion to the atom magnetic moment; 

• the positivity constraints; 

• the existence of an electron - positron sea and its role in the sum 
rules. 

As scaling variable we take the null-plane momentum of the electron 
measured in the atom rest frame, 

k'^ = {ko + kz)rest frame = {kz/ PA,z)inf. mom. frame = ^^j 

(1) 

We prefer it to the Bj0rken variable xbj which is very small and 
depends on the nucleus mass. The kinematical limit for is Matom 
but typical values are \k+ — me] ~ Zamg. 

We hope in this study to get a better insight of relativistic and 
spin effects in hadronic physics. The infinite momentum or null-plane 
descrition can also be interesting in atomic physics itself, since "deep 
inelastic" experiments can also be made with atoms, in particular: 

• Compton profile measurements: 

7(i^) -|- bound e~ — >■ ^{K + Q) + free e~{k'), 
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• Moeller or Bhabha scattering: 

e^{K) + bound e~ ^ e^{K + Q) + free e-{k'), 

• annihilation: e'^{K) + bound e~ ^{K + Q) + jik). 

The Mandelstam variables s = {K+Q + k'f, t = Q^ and u = {K-k'f 
are supposed to be large compared to rrig. A ten MeV beam is suffi- 
cient for that. We take the z axis opposite to the beam direction. In 
the laboratory frame the final particles are ultrarelativistic nearly in 
the — z direction. The components k'^ and k^^ of the electron momen- 
tum just before the collision are given by 

fc+ ~ Q+ , (2) 

ky ~ —P'rp(nudeus) = k^ + Qt- (3) 

k+ can be measured with one detector, k~^ and ky need two detectors. 
The definition ^ of ky is ambiguous due to the final state Coulomb 
interaction. 



2 Joint (b, k^) distribution 

Being observable quantities, the operators k^ and ky should be de- 
fined in the gauge independent way 

k+ = ido-V{r)-id,-A,{r), kr = -iVrir) - At- (4) 

They do not commute: [/c''',kT] = — iVT^(r), where — Vry(r) is the 
transverse part of the Coulomb force. Therefore one cannot define 
a joint distribution q{k~^ ,\i.T) in an unambiguous way. Leaving this 
problem for the next section, we can at least define the joint distribu- 
tion q{k~^ ,h) in the impact parameter representation. This quantity 
plays a role in double H + H collisions in which both nuclei and both 
electrons collide. 

From the known the Dirac wave function of the hydrogen atom [1] 
vl/(t,r) = ^(r)e-^^*, (5) 
we can define the two-component null-plane wave function in b and 

/ + 00 
dz exi>{-ik+z + iEz-ix{h,z)} $(r), (6) 
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sinh ^ ( —] — sinh ^ ( -r- 



X(b,2;) = / dz' V{x,y,z') = -Za 

(8) 

The "gauge link" exp{— ix(b, z)} transforms ^ in the Coulomb gauge 
to ^ in the null plane gauge = {A^ = in the infinite momentum 
frame). The choice of zq corresponds to a residual gauge freedom. The 
quantity 

will be temporarily interpreted as the electron distribution in the 
atom. One has indeed 

However a significant re-interpretation will be given in Section 5. 

The gauge link makes q{k~^ ,h) invariant under a gauge transfor- 
mation, for instance V{r) — > y(r)-|- Constant, E — > E+ Constant. 
Such a shift of the potential is practically realized when electrons are 
added in far outer shells. Intuitively, this addition does not change 
the momentum distribution of the deeply bound electrons. 



3 Joint (ky, k^) distribution 

Notwithstanding the non-commutativity mentioned earlier, one can 
make a transversal Fourier transform of ([6]), 

<l>(kT,k+) = J d^h e-'^^-^ $(b,A;+) (11) 

and define a longitudinal-transverse momentum distribution 

q{kT, k+) = cDt(kT, k+) cD(kT, k+), (12) 

normalized to 

q{k+)= Jq{kT,k+)d'-kT/{27rf. (13) 

q{kT,k~^) depends on zq, which is a remnent of the ambiguity. How- 
ever, an appropriate choice of zq turns this apparent disease into an 
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advantage [31 [H [5]. Taking zq = — oo for the Compton reaction, the 
factor e"*'^^'''^^ in ([6]) just takes care of the final state interaction: it 
describes the distortion of the scattered electron wave function by the 
Coulomb potential, in the eikonal approximation. Similarly, taking 
zq = +00 for the annihilation reaction, it describes the distortion of 
the initial positron wave function. Thus q(kT, which depends on 
Zq, has no precise intrinsic character. One can just consider a "most 
intrisic" definition with zq = 0. 



4 Spin dependence of the electron den- 
sity 



In formulas ()5l ll2p the angular momentum state of the atom was not 
specified. We assume that the electron is in the fundamental n = 1, 
j = 1/2 state and the nucleus is spinless. Let S"^ = 2(j) and S*^ = 2(s) 
denote the atom and electron polarization vectors. The unpolarized 
electron density in (b, space in a fully polarized (|S"^| = 1) atom 
is 

q{h,k+;S^) = $T(b,A;+;S^) $(b,A;+;S^) (14) 
and the electron polarisation is given by 

S^(b,A;+;S^) q{h,k+;S^) = ^^{h,k+;S^) a <^{h,k+;S^) . (15) 

Taking into account parity and angular momentum conservations, the 
density of electrons with polarization S*^ in a polarized atom can be 
written as 

(7(b,fc+,S^;S^) = ^i^^[l+Co„ S^-n+C„o S'-n+Cnn (S^-n)(S^-n) 

+ Cu StS^+Ci^ StiS'^-TT) +C^i iS^-TT)S^+C^^ (S^-7r)(S^-7r)], (16) 

where tt = b/6 and fi = z x vr. The Cj^'s also are functions of b 
and k'^ . A similar equation can be written in the representation. 
Integrating (jl6p over b leaves the following spin correlations: 

g(fc+,S^;S^) = ^ [q{k+) + Aq{k+)SIS^ + 6q{k+) • 

(17) 

where Ag(A;+) and Sqlk'^) are the helicity and transversity distribu- 
tions. 



5 



4.1 Formulas for the polarized densities in 

(b, k+) and (ky, k+) 

For the jz = +1/2 state, ^ of Eqs.([6l fTT]) can be written as 

$(b, k^; = +z) = _ J^,, ^ , $(kT, = +z) = _ 

(18) 

For the jz = — 1/2 state, 

c|>(b, = -z) = ( ) , $(k„ = -z) = ( 

(19) 

For other orientations of S"^, one takes Hnear combinations of (jlSp and 
(|19p . The (b, k~^) distribution depends of zq only in an over-all phase. 
Chosing zq = 0, v{h, k~^) and w{h, k'^) are real and given by 

{l)-iy{r + ) m/r. (20) 

where ^ = Za/{1 + 7), 7 = E/nie = \/l — (Za)'^ and 
/(O = ' (2meZa)^+i/2 r^-i exp(-meZar) (21) 

is the IS radial wave function. Then, 

q{b,k^) = -w^+'y^ 

C„„(6,A;+) = 1 

Conib,k~^) = Cno{b,k+) = -2wu/iw^ + v^) 

Cu{b,k-^) = C^^{b,k+) = {w^-v^)/(w^+v^) 

CUb,k+)=C^i{b,k+) = 0. (22) 

Note that CQn{b,k^) ^ 0, which gives an asymmetrical impact pa- 
rameter profile for a transversely polarized atom. 

The {hT,k~^) distribution depends on zq. Taking zq = ^00 makes 
([5D diverges. In practice we will assume that l^ol is large but finite, 
accounting for a screening of the Coulomb potential. It gives 

X(b, z) = -Za \ smh-\z/b) - e{zo) ln(2|zo|/6) 1 , (23) 
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with e(0) = and e(=Foo) = ^1, the upper sign corresponding to 
Compton scattering and the lower sign to annihilation. Modulo an 
overall phase, 

* U 2vr /~ 6 d6 6'^-(^o) ( Mkrb) w{b k+) \ 
v) Jo \ Ji{kTb) v{b,k+) J 



The analogue of (122]) is 

q{kT,k~^) = \w\'^ + \v\'^ 

Cnn{kT,k+) = 1 
Con{kT,k+)=Cno{kT,k+) = 2^{v* w) / {\w\^ + 
Cu{kT,k+)=CUkT,k-^) = (|^|2_|i)|2)/([^[2 + |^[2) 

CUkT,k-^) = -C^i{kT,k~^) = 2^{v*w)/{\w\^ + \v\^). (25) 

These coefficients are related to the structure functions listed in Ref . [6] . 
For the "most intrisic" gauge zq = 0, w and v are real so that 
Conikr, k^) = (no Sivers effects). This is in accordance with time 
reversal invariance [5]. For the "Compton" and "annihilation" gauges 
(zo = =f1); w and v are complex numbers, so that Sivers [7j effect 
lCon{kT,k+) / and the Boer-Mulders [6] effect (C„o(fer,A:+) / 0) 
take place. 

In the Compton case the factor b~^^°' behaves like a converging 
cylindrical wave. Multiplying <&(r), it operates as a boost toward the z 
axis, interpreted as the "focusing" of the final particle by the Coulomb 
field [1] . This focusing converts the asymmetry in b for a transversely 
polarized atom into the Sivers asymmetry in kj-. The opposite effect 
(defocusing of the positron) takes place in the annihilation case. 



4.2 Sum rules 

Integrating ()17p over k^ , one obtains the vector, axial and tensor 
charges 

q = j ¥{r- S^) ^(r; S^) = 1, (26) 
Aq = S^-J d^v vl/t(r; S^) S ^(r; S^) = ^Y^' (^7) 

5q = S^- j d\ ¥{r- S^) (i S ^(r; S^) = ^^^p^- (28) 

Note the big "helicity crisis", Ag = 1/3 instead of 1 as naively ex- 
pected, for Za = 1. 
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4.2.1 Sum rule for the atom magnetic moment 

Consider a classical object of mass M, charge Q, spin J and time- 
averaged magnetic moment /2 in its rest frame. In this frame, the 
centre of energy re and the average center of charge (rc) coincide, 
say at r = 0. Upon a boost of velocity v, rc and (rc) undergo the 
lateral displacements 

be = V X J/M , (be) = V X . (29) 

he and (be) coincide if the gyromagnetic ratio has the Dirac value 
Q/M. In our case, be is negligible due to the large nucleus mass, 
therefore the magnetic moment is almost totally anomalous. In the 
infinite momentum or null-plane frame (v ~ z) one observes an electric 
dipole moment jSj 

- e (b) = fiatom Z X S^, (30) 

which we can calculate from Con{b, k^). Weighting (fT4|) with bx = x 
for S"^ = y one obtains 

(x) = -(l + 27)/(6me) (31) 

which is accordance with the relativistic result for the atomic magnetic 
moment ha = — e (1 + 2'y) / (Qnie) (ignoring the anomalous magnetic 
moment of the electron itself). 

4.3 Positivity constraints 

The spin correlations between the electron and the atom can be en- 
coded in a positive-definite "grand density matrix" R [9], 

R = C^u ® Kf ■ (32) 

Here /u and v run from to 3, summation is understood over repeated 
indices, = I and Cqo = 1. R can be seen as the density matrix of 
the final state in the crossed reaction nucleus atom(S'^) + e^(— S*^). 
Besides the trivial conditions \Cij\ < 1 the positivity of R gives 

(1 ± Cnn? > {CnO ± Conf + {Ql ± C^^f + {C^l T Ci^f ■ (33) 

These two inequalities agree with those of Ref.[10]. Together with 
\Cii\ < 1 they are saturated by (f22l) or (f25l) . This maximal strength of 
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the spin correlation means that the information contained in the atom 
polarization is fully transferred to the electron, once the other degrees 
of freedom (k~^ and b or k^) have been fixed. If we integrate over 
k^, for instance, some information is lost and some positivity condi- 
tions get non-saturated. The same happens if there are "spectators" 
electrons which keep part of the information for themselves. 

After integration over b or kj^, we are left with the Soffer inequality 

2\6q{k+)\<q{k+) + Aq{k+), (34) 

which are saturated by (12611280 . 

Note that a complete anti-correlation between the atom and the 
electron spins, Cu = Cnn = C'tttt = — 1 and Cj^j = 0, leading to 

g(b,A;+,S^;S^) = g(b,A;+) (1 - • S^)/2, 

violates the positivity conditions, although the last expression is pos- 
itive for any and S"^. In fact such a correlation would make 
{A, e'^\R\A, e^) negative for some entangled sisXes \ A, e"*") in the crossed 
channel, in particular the spin-singlet state [9]. 

5 The electron-positron sea 

The charge rule (|26p receives positive contributions from both posi- 
tive and negative values of k~^ . So the contribution of the positive A;"*" 
domain is less than unity. On the other hand, physical electrons have 
positive k'^ . It seems therefore that there is less than one physical elec- 
tron in the atom. This paradox is solved by the second quantification 
and the introduction of the electron-positron sea. 

Let us denote by |n) an electron state in the Coulomb field. Quan- 
tizing the states in a box, we take n to be integer. Negative n's are as- 
signed to negative energy states. Positive n's up to ns label the bound 
states {—me < En < +We) and the remaining ones from ub + 1 to 
-|-oo label the positive energy scattering states. En > +771^. Let \k,s) 
be the plane wave with four-momentum k and spin s, solution of the 
free Dirac equation. The destruction and creation operators in the 
interacting and free bases are related by 

ak,.s = {k,s\n) an, = ^ aj, (n|A;, s). (35) 

n n 
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In the Dirac hole theory, the hydrogen-hke atom is in the Fock state 
\Hn) = On '^-2 ■ ■ ■ ^-oo |Dirac-bare nucleus). (36) 

"Dirac-bare" means that all Dirac states, including the negative en- 
ergy ones, are empty. The number of electrons of momentum k and 
spin s (with positive and k'^) in the atom is 

Ktom(k,s) = (4,5 OLk,s)atom = \{k,s\n)\'^ + ^ \{k,s\n')\^. (37) 

n'<0 

A stripped ion is a "Dirac-dressed" nucleus, all negative energy states 
being occupied. For the ion the factor of p6|) is missing and the 
first term of (j37p is absent. By difference. 



\{k,s\n)\'= f ^q{k+), (38) 
7Tr„ . Jk+>0 2,11 



fc,fcO>0 s 



the last expression being for the continuum limit (A;, s\n) — > <I?(kr, fc"*"). 

Positrons are holes in the Dirac sea. The number of positrons of 
momentum k (with positive k^ and k'^) and spin s in the atom is 

Ktom(k,s) = (a_fc^„ga^__)atom= E \{-k,-s\n')\\ (39) 
For the ion, the condition n' 7^ n is relaxed. By difference, 

NL-Ktm= -^Qik^), (40) 

Jk+<o ^vr 

where we have made the change of variable k = —k. The sum rule 
(|26p can therefore be interpreted in the following way: 

• for k^ > 0, q{k^) = (e~ distrib. in atom) - (e~ distrib. in ion) 

• for k^ < 0, qik'^) = (e^ distrib. in ion) - (e"^ distrib. in atom) 
Thus 

[ — N'^~ \ -\- ( N^'^ — N^^ ^ = 1 (A^\ 

\-^^atom ^^lonj ^ l^'ion atom J V^^l 

where each braket € [0, 1]. Introducing Q'^ = —N^^, this can be 
rewritten as 

Qatom ~ Qion ~ ^ (^2) 

Qion is the electronic charge renormalisation of the ion on the null 
plane. It is more likely positive, maybe infinite for a pointlike nucleus. 
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The renormalisation Qatom — 1 of the atom is equal to it. It may be 
interesting to relate Qfon with the result of covariant QED. 

dKLm/dk+ , dN±Jdk+ , dNf-Jdk+ and dN^Jdh^ are sepa- 
rately measured in the deep inelastic reactions listed in Section 1 and 
their generalizations to the electron - positron sea, for instance 

-i{K) + sea ^ -i{K + Q) + free e^{k'). (43) 

A "sea" electron can be equally understood in the sense given by 
Feyman in the parton model or by Dirac in the hole theory. It gives 
the second term of ()37p and the whole process is 

J + A ^ -f' + A + slow e'^ + fast e~ . (44) 

A "sea" positron is understood in the Feynman sense only. In the hole 
theory, an electron of large negative energy is lifted to a bound or slow 
free state |n'). It gives the right-hand side of (j39p . The whole process 
is 

j + Hn^ + fast e+ + y. (45) 

Of course, one can permute the roles of the electron and the positron 
in the Dirac theory; then Feynman and Dirac sea positrons become 
equivalent. 



6 Momentum sum rule 

q{k+) obeys a momentum sum rule which, like the charge sum rule, 
applies to the difference between the atom and the ion. 

The null-plane momentum of the ion (=nucleus) can be decom- 
posed into a matter part and a Coulomb field part: 

= M^ = P+^,,^^ ^ + {E^} . (46) 

^matter N includes the momentum of the electron cloud which renor- 
mahze the ion charge. P^-^^^ {Eat} is the flux of the T^'" = T^" + T^'" 
component of the energy-momentum tensor T^^ {Ejv} of the nucleus 
Coulomb field E7v(r) across the null plane: 

PfieiA^N} = j {t'^{^n} + T^''{^n}) da. (47) 

= jdxdydz (t°° + + + T^^) (48) 

= j dx dy dz (eI + (49) 
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We have used duy = (1,0,0,1) dx dy dz, = ^(E^ + B^), T^^ = 
^00 _ £;2 _ j^oz ^ j^zo ^ ^^^^ _ ^^^^^ Similarly, for the atom, 

we have 

= Ma = Patter N + PLe e + + ; Be} . (50) 

the electron magnetic field being included. Here P^^^^, g, Eg and Bg 
take only into account the difference between the atomic and ionic 
electron clouds. Substracting ([^61) from ([50]) . 

E^Ma-Mn = P^ +Ptnf (51) 

P^f- results from the crossed terms in Eat and Eg or Bg of T^'^ . Its 
value is 

Pm ~ j [ ^^^^ ~^ ^'^y^ ~^ ^-^^y ~ ^'^^^ 1 

= ^y'd3rEjvEg = ^(y). (52) 

(y) is the average potential energy. The terms in E/^rBg have disap- 
peared upon angular integration. 

Pg"*" = -P^^^g g + P^g;^{Eg} is the mean value of the null-plane 
mechanical momentum /c"*" of the physical electron (more precisely the 
"atom - minus - ion" part of it). Inserting = E — idz in Eqs.®, 
one obtains 

_^ ^ q{k+) = E - ^{V) (53) 

with 

(V) = J d^r ^t(r) V{r) ^(r) = -me{Zaf/j. (54) 
Eqs. (l51ll53l) constitute the momentum sum rule. 

7 CONCLUSION 

This study has shown the rich spin and structure of the hydrogen- 
like atom at large Z when it is observed in the infinite momentum (or 
null-plane) frame. Without the complications of QCD, like gluon self- 
interaction and confinement, many properties attributed to the leading 
twist hadronic structure functions have been found and clearly inter- 
preted here, in particular: the sum rules, the spin crisis, the connection 
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between (b) 7^ and the Sivers effect, the relation between (b) and the 
magnetic moment, the role of spectators in the positivity constraints, 
the existence of a Feynman sea. With this "theoretical laboratory" 
one may also investigate non-leading twist structure functions, elastic 
form factor a la Isgur-Wise, etc. Our results are interesting also in 
pure QED. We have seen a connection between the nucleus charge 
renormalization and the unpolarized deep inelastic structure function 
of the cloud of a stripped ion target. Thus the charge renormal- 
ization can be analyzed experimentally in deep inelastic Compton, 
M0ller or annihilation processes. The same reactions at the 10 MeV 
energy scale can test the relativistic corrections to the electronic wave 
functions of large Z atoms. 

The powerpoint document presented at QEDSP-2006, which con- 
tains figures not presented here, can be obtained upon request to one 
author (X. A.). 
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